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Recent work has shown that differential rotation, producing large scale drifts of fluid elements 
along stellar latitudes, is an unavoidable feature of r-modes in the nonlinear theory. We investigate 
the role of this differential rotation in the evolution of the I = 2 r-mode instability of a newly 
born, hot, rapidly rotating neutron star. It is shown that the amplitude of the r-mode saturates a 
few hundred seconds after the mode instability sets in. The saturation amplitude depends on the 
amount of differential rotation at the time the instability becomes active and can take values much 
smaller than unity. It is also shown that, independently of the saturation amplitude of the mode, 
the star spins down to rotation rates that are comparable to the inferred initial rotation rates of the 
fastest pulsars associated with supernova remnants. Finally, it is shown that, when the drift of fluid 
elements at the time the instability sets in is significant, most of the initial angular momentum of 
the star is transferred to the r-mode and, consequently, almost none is carried away by gravitational 
radiation. 

PACS numbers: 04.40.Dg, 95.30.Lz, 97.10.Sj, 97.10.Kc 



I. INTRODUCTION 

The interest in r-modes, first studied more than twenty 
years agojlj , has increased dramatically when it was dis- 
covered U/and afterwards confirmed more generally Q , 
that r-modes are driven unstable by gravitational radi- 
ation reaction in perfect-fluid stars with arbitrary small 
rotation. Soon afterwards, it was shown 0, that in a 
newly born, hot, rapidly rotating neutron star, bulk and 
shear viscosity do not suppress the r-mode instability for 
a wide range of relevant temperatures and angular veloc- 
ities of the star. As a result, the neutron star could spin 
down to just a small fraction of its initial angular velocity, 
thus providing a possible explanation for the relatively 
small spin rates of young pulsars in supernova remnants 
0,0 . The gravitational waves emitted by the young neu- 
tron star during this spin-down phase could be detected 
by enhanced versions of laser interferometer detectors Q . 
It was also suggested that in accreting neutron stars in 
low-mass x-ray binaries the gravitational-wave emission 
due to the r-mode instability could balance the spin-up 
torque due to accretion, thus limiting the maximum an- 
gular velocity of these stars to values consistent with ob- 
servations 

Of fundamental importance in judging the astrophysi- 
cal relevance of the r-mode instability is the determina- 
tion of the saturation amplitude of the mode. In recent 
years, several authors have addressed this issue, both an- 
alytically and numerically, taking into account different 
nonlinear effects HHElHIliElHm. 

An approximate analytical expression for differential 
rotation induced by r-modes was first derived using the 
linearized fluid equations by expanding the velocity of 
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a fluid clement located at a certain point in powers of 
the mode's amplitude, averaging over a gyration, and re- 
taining only the lowest-order nonvanishing term 0, llf)j | . 
This differential rotation can then interact with the mag- 
netic field of a newly born, hot, rapidly rotating neutron 
star, limiting the growth of the r-mode instability or, for 
strong magnetic fields, even preventing it from develop- 
mg HEl 

The first numerical study of nonlinear r-modes was 
performed for a rapidly-rotatin g re lativistic star without 
a gravitational radiation force [lj- These studies show 
there is no suppression of the mode, even when its initial 
amplitude is of order unity, and also confirm the exis- 
tence of differential rotation induced by r-modes. In a 
subsequent numerical study of the nonlinear evolution 
of r-modes of a rapidly-rotating Newtonian star with a 
gravitational radiation force, it was shown that the am- 
plitude of the mode grows to order unity before strong 
shocks near the stellar surface rapidly damp the mode 
Il4j. However, a more recent numerical simulation 
has found no evidence that the decay of the mode's am- 
plitude is due to such shocks near the surface of the star; 
instead, the catastrophic decay of a r-mode's amplitude 
of order unity is due to a leaking of energy into other fluid 
modes, leading to a differentially rotating configuration 

Recently, the nonlinear coupling between stellar iner- 
tial modes has been analyzed, with the conclusion that 
r-modes may saturate at much lower values than previ- 
ous investigations had revealed fl|j. Despite this fact, it 
was found that the r-mode instability could still explain 
the spin clustering at the fast end of the spin distribu- 
tion of neutron stars in low-mass x-ray binaries and that 
gravitational waves from newly born, hot, rapidly rotat- 
ing neutron stars, as well as from old neutron stars in 
low-mass x-ray binaries, could be detected by enhanced 
versions of laser interferometer gravitational wave detec- 
tors 0. 
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In this paper, we are concerned with the role of differ- 
ential rotation in the evolution of the r-mode instability. 

Recently, a nonlinear extension of the linear r-mode 
perturbation was found within the nonlinear theory up 
to second order in the mode's amplitude a in the case 
of a Newtonian, barotropic, perfect-fluid star rotating 
with constant angular velocity O p"?j | . This nonlinear ex- 
tension describes differential rotation of pure kinematic 
nature that produces large scale drifts along stellar lat- 
itudes. This solution contains two separate pieces, one 
induced by first-order quantities and another determined 
by the choice of initial data. These two pieces cannot 
cancel each other, since one is stratified on cylinders and 
the other not. Thus, differential rotation is an unavoid- 
able kinematic feature of r-modes. 

The differential rotation found in Ref. [jjj does con- 
tribute to the second-order physical angular momen- 
tum of the r-mode, implying that the phenomenological 
model proposed in Ref. |6j to study the evolution of the 
r-modc instability in newly born, hot, rapidly rotating 
neutron stars has to be modified so that the physical an- 
gular momentum of the perturbation includes, not only 
the canonical angular momentum, but also the contri- 
bution from second-order differential rotation. In this 
paper, we study such a modified model, arriving at the 
conclusion that differential rotation plays an important 
role in the evolution of the r-mode instability. 

In section^we derive an analytical expression for the 
physical angular momentum of the r-mode perturbation, 
using the second-order solution of Ref. [l7|. In section 
IIHI the model for the evolution of the r-mode instability 
is presented and a system of nonlinear, coupled, differen- 
tial equations for the mode's amplitude a and the star's 
angular velocity is derived. In section llVl the analyti- 
cal solution of the above mentioned system of differential 
equations is presented, revealing that differential rota- 
tion does play an important role in the evolution of the 
r-mode instability. Finally, in section we present the 
conclusions. 



II. PHYSICAL ANGULAR MOMENTUM OF 
THE .R-MODE PERTURBATION 

For small-amplitude perturbations of slowly rotating, 
Newtonian, barotropic, perfect-fluid stars, the r-mode 
solution to the linearized Euler and continuity equations 
is, in spherical coordinates (r, 6, (f>), given by: 

5^v r = 0, (la) 
8 {1) v e = aQQl(^j 1 sin 1 - 1 8sm(l<j) + ujt), (lb) 

<5 (1 V = a&Ctli^) sm'- 2 0cos0cos(^ + u;t), (lc) 

where 8^v l are the contravariant components of the 
first-order Eulerian change in velocity, R and f2 are, re- 
spectively, the radius and the angular velocity of the un- 



perturbed star, 

u = -m + — (2) 

is the mode's angular frequency in an inertial frame, a is 
the mode's amplitude and C/ is chosen to be 

To the velocity field JQ) corresponds a Lagrangian vector 
displacement 

C (1)r = 0, (4a) 

i {1)e = --aCil(l + l) (-Y 1 sin 1 - 1 9cos{l<j) + Ljt), 
2 \ R/ 

(4b) 

^ = l a Ql(l + 1) 1 sin 1 ' 2 6 cos9 sin(l<j) + Lut), 

2 V RJ 

(4c) 

satisfying, at lowest order in J7, the equations 

5^v i = d t £, {1)l + v k V k ^ 1)l - i [1)k V k v L (5) 

and 

V^W*) - 0, (6) 

where v l = fWJ, and p are, respectively, the fluid veloc- 
ity and the mass density of the unperturbed star. The 
knowledge of the first-order vector displacement is 
sufficient to compute the canonical angular momentum of 
the r-mode perturbation [18|; using Eq. |0J, one obtains: 

J c = - j P d^(d t £Uv k v k ^)dv 

= -^a 2 ni(l + l)R 2 - 21 j\r 2l+2 dr. (7) 

However, the canonical angular momentum J c is not the 
full physical angular momentum at second order d^ 2 'J, 
the difference being given by 0] 

5 {2) J- J c = ^Jpv i Af ) v i dV, (8) 

where A^Vi is the second-order Lagrangian change in 
velocity 

+ ^ (2) + , fc (v4 2) +V^ 2) ), (9) 

which, as opposed to the canonical angular momentum, 
also contains terms linear in the second-order Lagrangian 
displacement vector ^ 2 ^ a , in addition to terms quadratic 
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in ^ 1 " >a . Therefore, the computation of the physical an- 
gular momentum at second order requires the knowledge 

Of £(2)0. 

In a recent paper [l7j . a nonlinear extension of the 
linear r-mode was found within the nonlinear theory up 
to second order in the mode's amplitude in the case of 
a slowly rotating, Newtonian, barotropic, perfect-fluid 
star: 



S^v r = 0, 
S^v e = 0, 



1 



r \ 



21-2 



-cfnCW-m-) sin 
+ a 2 nAr N - 1 sm N - 1 e, 



21-4 , 



(10a) 
(10b) 

9 

(10c) 



where A and N are arbitrary constants fixed by initial 
data. This second-order solution represents differential 
rotation, producing large scale drifts of fluid elements 
along stellar latitudes. The first term on the right-hand 
side of Eq. i|10cfl is induced by first-order quantities, while 
the second term is determined by the choice of initial 
data. Since these two terms cannot cancel each other, 
differential rotation is an unavoidable feature of r-modes. 
Let us in the present work restrict ourselves to the case 
N = 21 — 1. We also redefine A to be 



A = -C?l 2 {l + l)R 
2 



2-21 



K. 



(11) 



where if is a constant fixed by the choice of initial data. 

At lowest order in fi, the contravariant components of 
the second-order Lagrangian displacement vector 
corresponding to the axisymmetric time-independent ve- 
locity field (|10|) . are determined by the equations 01 

- ^V^V, (12) 



(2)k 



1 



V fe £ (1)m V m £ 



(i)k 



C (2)fe v fe p = -^ (1)fc £ (1)m v fe v mP , 



(13) 
(14) 



which yield the solution [T^ 

= l_ a 2 C 2 l 2 [l + l) 2 R ^ 



/ r \ 21-1 



x sin 2 ^ 2 6» (sin 2 6» - 2) 



(15a) 



,{2)9 
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x sin 2/ ~ 3 9 cos 9 (sin 2 + 21 - 2) , (15b) 



x sin 



21-2 



9 t + C(r, ( 



(15c) 



where C is an arbitrary function of r and 9. In Eq. Ijl2|l - 
(|14fl . terms quadratic in first-order quantities give rise 
to double "frequency" terms of the type cos[2(l<fi + u>t)] 



inducing a second-order solution corresponding to an os- 
cillating response at a "frequency" twice that of the linear 
r-mode; however, in this article these double frequency 
oscillations will not be considered since they give a zero 
contribution to the second-order physical angular mo- 
mentum of the r-mode perturbation. Using the above ex- 
pressions for the second-order Lagrangian displacement 
vector £( 2 )' , the physical angular momentum at second 
order is computed to be 



8&J 



-a 2 n [2KI + 1)(2Z + 1)] R 2 - 21 



x / pr 
lo 



21+2 



dr. 



(16) 



As can be seen from Eqs. (0) and l|16|l . the physical 
angular momentum of the r-mode perturbation is not 
equal, in general, to the canonical angular momentum. 
The former contains also a part linear in the second- 
order Lagrangian change in velocity, which, as pointed 
out in Ref. 18], is related to conservation of circulation 
in the fluid. Since, as shown in Ref. [TtI ]. at second or- 
der r-modes do not conserve vorticity, it follows that, in 
general, the physical and canonical angular momentum 
do not coincide. However, specific choices of K can be 
made such that the integral in Eq. JSJ vanishes. Such a 
case (K = —2, for I = 2) was studied in Ref. 0] within a 
phenomenological model for the evolution of the r-mode 
instability. There is not, however, to our knowledge, a 
physical condition that forces K to take such a particular 
value. Therefore, in this paper we study the evolution of 
the r-mode instability for arbitrary values of K. 



III. THE EVOLUTION MODEL 

In this paper, the simple model proposed in Ref. 6] to 
study the evolution of the r-mode perturbation in newly 
born, hot, rapidly rotating neutron stars is adopted with 
two modifications. First, the physical angular momen- 
tum of the r-mode perturbation is not taken to be just 
the canonical angular momentum as in Ref. 0- Instead, 
we will use the full physical angular momentum l|16|l . 
which, as already discussed in the previous section, be- 
cause of the presence of differential rotation induced by 
r-modes, is different, in general, from the canonical angu- 
lar momentum. A second modification, less significant, 
is related to the proposal of Ref. to deduce the evolu- 
tion equations for the mode's amplitude a and the star's 
angular velocity f2 just from angular momentum consid- 
erations. 

Since the most unstable r-mode is the I = 2 mode, we 
will focus in this paper our attention only on the evolu- 
tion of this mode. We will also assume, as in Ref. 0, 
that the mass density p and the pressure p of the fluid 
are related by a polytropic equation of state p = kp 2 , 
with k such that M = IAMq and R = 12.53 km. 
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In our model, it is assumed that the total angular mo- 
mentum of the star is given by 



J = m + 8 {2) J 



(17) 



where the momentum of inertia of the equilibrium con- 
figuration I is given by 



I = 



8tt 



pr 4 dr = I MR 2 , 



(18) 



with I = 0.261, and the physical angular momentum of 
the r-mode perturbation J is given by 



1 i r R 

-a 2 n(4K + 5)—2 / pr 6 dr 

2 R Jo 



-a 2 n(4:K + 5)JMR 2 , 



(19) 



with J = 1.635 x 1(T 2 . 

The total angular momentum of the star J = J(a, Q) 
decreases due to the emission of gravitational radiation 
according to |(| 



^=3JMR 2 ^. 
at t G r 



(20) 



where the gravitational-radiation timescale tgr is given 
by i 



1 

TGR 



1 

TGR 



n 2 

irGp 



(21) 



with the fiducial timescale tgr = —3.26 s. In the above 
equation, p is the average mass density of the star and G 
is the gravitational constant. 

From Eqs. i|17l) (1211) it is straightforward to obtain a 
differential equation for the time evolution of a and f2: 



1 



l + -(4K + 5)Qef 



~dt 



] -(4K + 5)Qna^ = 2Q— 
> at tgr 



(22) 



where Q = 3 J/ (27) = 0.094. 

Following the proposal of Ref. to deduce the evo- 
lution equations just from angular momentum considera- 
tions, we assume that the physical angular momentum of 
the r-mode perturbation 5^ J increases due to the emis- 
sion of gravitational radiation and decreases due to the 
dissipative effect of viscosity, 



dd^J 
dt 



-26® J ( — 

\ TGR 



1 

TV 



(23) 



For the viscous timescale ry, we take the expression de- 
rived in Ref. |4( for the simple case of the linear r-mode 
of a newly born, hot, rapidly rotating neutron star 
with shear and bulk viscosity: 



1 

TV 



1 fl0 9 K 

fs 



T 



n 2 



-(— V, 

f B V10 9 K7 \irGp 



(24) 



with the fiducial timescales fs — 2.52 x 10 s s and tr — 
6.99 x 10 s s. Several authors have taken into account 
other dissipative mechanisms, but in this paper we shall 
restrict ourselves, for illustrative purposes, to the above 
expression. 

From Eq. 123fl one obtains, using Eq. I|19|) . a second 
differential equation for the time evolution of a and CI: 



„_da dQ, 
2n—+a — 
dt dt 



\TGR 



1 

TV 



(25) 



Note that in Ref. 6] it was assumed that it was the energy 
of the perturbation (in the rotating frame) , and not the 
angular momentum, that increases due to emission of 
gravitational waves and decreases due to viscosity. For a 
constant these two approaches are coincident; however, 
for a varying they differ in quantities of the order of 
Qa 2 |Hj. 

From Eqs. 11221) and ll'ol) it is straightforward to obtain 
a system of two, first-order, coupled, differential equa- 
tions determining the time evolution of the amplitude of 
the r-mode a(t) and of the angular velocity of the star 

n(t)-. 



dn 

~dt 
da 
~dt 



3 TGR 



:W + 5)Q^! 

3 tv 



l + -(K + 2)Qa 2 



1 + -(AK + 5)Qa 2 



a 
tgr 
a 

TV ' 



(26) 



(27) 



For K = —2, the above equations coincide with 
Eqs. (2.7) and (2.8) of Ref. [Hf (for l/r M = 0) and 
with Eqs. (3.14) and (3.15) of Ref. (up to terms of 
order Qa 2 ). In this paper, as already mentioned above, 
we study the evolution of the r-mode instability for an 
arbitrary value of K. 



IV. THE ROLE OF DIFFERENTIAL ROTATION 
IN THE EVOLUTION OF iJ-MODES 

The condition Tqr{CI) + r v 7 1 (0,T) = gives the sta- 
bility curve, i.e., the set of points in a diagram (f2, T) for 
which the damping effect of viscosity balances exactly 
the driving effect of gravitational radiation. For a newly 
born, hot, rapidly rotating neutron star, there is an in- 
terval of temperatures and angular velocities of the star 
for which the gravitational timescale tqr is much smaller 
than the viscous timescale tv, implying that the evolu- 
tion of a and 17 is determined, in a good approximation, 
by the equations: 



dfl 
~dt 
da 
~dt 



1 



2)Q- 



tgr 



-{K + 2)Qa 2 



a 

TGR 



(28) 
(29) 



5 



These equations remain a good approximation to 
Eqs. and (|77|) for a period of time during which the 
temperature and angular velocity of the neutron star are 
such that the corresponding point, in a (f2, T) diagram, 
lies well above the stability curve. For the model we have 
been considering, this period of time is about one year 
□ 

In solving the above system of equations, initial condi- 
tions at t = to are chosen such that |<5^ 2 ^ J(to)\ <C Ifl(t ). 
For a = a(t Q ) = 1(T 6 , this implies that |A| « 10 13 . 
In the case A < —5/4, for which <r 2 )j < 0, as the 
amplitude of the mode grows due to the gravitational- 
radiation instability, the total angular momentum of the 
star, given by Eq. (|17J) . decreases and eventually becomes 
negative. To avoid this, the growth of the mode's ampli- 
tude has to be stopped by hand at a saturation value 
a sa t ^ \J— 3/[(4A + 5)Q]; integration is then carried on 
with a new set of equations for which a — a sa t and the 
evolution of the angular velocity Q is determined from 
Eq. with da/dt = 0, 



dfl 
~dt 



2Q 



T GR l + l(4K + 5)a* at ( 



(30) 



In the case A ^ —5/4, for which 5^ J ^ 0, as the ampli- 
tude of the mode grows, the total angular momentum de- 
creases but remains always positive and, therefore, there 
is no need to saturate the mode by hand. In this case, 
a(t) and Q(t) are determined solely by Eqs. (|2*5|) and i(2*§|) 
and, as will be shown below, the mode saturates naturally 
at a value that depends on the parameter A, namely, 
a sat oc (A + 2)~ 1 / 2 . As already mentioned above, the 
parameter A, introduced in Eq. is fixed by initial 

data, giving the initial amount of differential rotation as- 
sociated to the r-mode. Thus, depending on the choice 
of initial data, the saturation amplitude of the r-mode 
can be of order unity, as assumed in the first models of 
evolution of the r-mode instability Q , or very small, as 
recent studies within the nonlinear theory seem to imply 

M 

The requirement that the total angular momentum is 
always non-negative, implies that the right-hand side of 
Eq. H29f) is always positive, i.e., a(t) increases in time 
for any value of A. It is also worth mentioning that, 
according to Eq. (|28[l . for A ^ —2, f2(i) evolves on a 
gravitational timescale, that is, a change in the back- 
ground motion of the star occurs simultaneously to the 
growth of the mode due to the gravitational-wave insta- 
bility. This result is not completely unexpected, since 
such kind of behavior has already been reported in the 
case of a toy model of a thin spherical shell of a rotating 
incompressible fluid [20|. 

Let us first consider the case A ^ —5/4. As already 
mentioned above, in this case the total angular momen- 
tum of the star is always positive and the amplitude of 
the mode a(t) and the angular velocity of the star f2(t) 
are determined solely by the system of Eqs. I|28|) and (12911 . 



from which it is straightforward to obtain: 
dfl f(A + 2)Qa 



-da. 



n l + f(A + 2)Qa 2 
This differential equation yields the solution 
l + |(A + 2)Qa 2 



l + |(A + 2)Qa 2 ' 



(31) 



(32) 



where fio = Q(to) and ao = a(to) are, respectively, the 
initial angular velocity of the star and the initial ampli- 
tude of the r-mode perturbation. Using the above ex- 
pression to eliminate f2 from Eq. (|29|l and integrating, 
one obtains: 



1 

TGR 



n 



l + -(A + 2)Qa 2 



(t-t ) 



In- 



5! 



Q'O 



x a 



^ 2n(-> - //)!//! 

n— 1 ^ 1 



; {K + 2)Q 



(33) 



In the initial stages of the evolution of the r-mode 
instability, the right-hand side of Eq. (|33f) is dominated 
by the first term and a increases exponentially, 



a{t) ~ a exp < 0.027 



n 



K 



(t - t a ) 



(34) 



as expected; for later times, the right-hand side of 
Eq. (|33|) is dominated by the last term (n = 5) and a 
increases very slowly as 



Q 



a(t) ~ 2.48 ( ^ ) (t - t ) 



3/5 



1/10 _ 



1 



Va 



(35) 



In Eqs . Q and (jSB) . t - tp is given in seconds, CIk = 
(2/3)^irGp is the Keplerian angular velocity at which the 
star starts shedding mass at the equator, we have used 
for the gravitational timescale the value tgr = —3.26 s, 
and we have assumed that A and ao are such that 
4(A + 2)Qag/3 <C 1. The smooth transition between 
the regimes l|34f) and 135fl occurs for t — to — few x 10 2 
seconds. 

The equation describes the time evolution of 

the amplitude of the mode between the moment the 
gravitational- wave instability sets in and the moment vis- 
cosity effects become dominant and start damping the 
mode. It reveals that saturation of the r-mode occurs 
just a few hundred seconds after the beginning of the 
exponential growth of a. Thus, within the model of evo- 
lution we are studying, which includes the nonlinear ef- 
fect of the differential rotation induced by r-modes, the 
mode's amplitude saturates in a natural way (see Fig.^l. 

The saturation value of the amplitude of the mode de- 
pends crucially on the parameter A, a sa t oc (A + 2)~ 1//2 . 
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FIG. 1: Time evolution of the amplitude of the r-mode a for 
different values of K. The saturation value of the amplitude 
of the mode depends crucially on the parameter K, which is 
related to the amount of initial differential rotation associated 
with the r-mode. The initial values of the amplitude of the 
mode and of the angular velocity of the star are, respectively, 
ao = 10~ 6 and flo = CIk- 



The initial amount of differential rotation associated with 
the r-mode can be minimized by choosing K appropri- 
ately, namely, K ~ 0, with the consequence that the am- 
plitude of the mode saturates at values of order unity. In 
this case, other nonlinear effects as, for instance, mode- 
mode interactions, that saturate r-modes at values much 
smaller than unity |l6j| . are more relevant than differen- 
tial rotation in what concerns the saturation amplitude 
of the mode. However, if the initial differential rotation 
associated to r-modes is significant, when compared with 
Q!q , then the saturation amplitude a sa t can be as small 
as 10~ 3 — 10~ 4 . In this case, differential rotation plays 
an important role in the saturation of r-modes. 

Let us now turn our attention to the time evolution of 
the angular velocity of the star fi. Integrating Eq. I|28jl. 
after eliminating a with Eq. (|32[l . one obtains the follow- 
ing solution: 



2 

TGR 



[1 



l + -(K + 2)Qal 



(t-to) 



5 

£ 

71=1 



l(K + 2)Qal] 



In 



n 



In- 



Oo 

n [ \ f] 

1 + §(^ + 2)Qq2- 
l(K + 2)Qa 2 



- 1 



(36) 



In the initial stages of the evolution of the r-mode 
instability, the right-hand side of Eq. (|36fl is dominated 
by the last term and f2 decreases as 



^)^l-l(^ + 2)Q^exp. 



0.054 



(37) 

for later times, the right-hand side of Eq. I|36|) is domi- 
nated by the first term (n — 5) and fl decreases slowly 



1.30 



n K 



-6/5 



(f - h) 



-1/5 



(38) 



In Eqs. (|37|l and l)38[l. again, t — to is given in seconds, we 
have used for the gravitational timescale the value tqr = 
—3.26 s, and we have assumed that K and ao are such 
that A(K-\-2)Qa 2 l /'& <C 1. The smooth transition between 
the regimes l|37|) and l|38|) occurs for t — to — few x 10 2 
seconds (see Fig- - Remarkably, in the later phase of 
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FIG. 2: Time evolution of the angular velocity of the star 
Q, for different values of K. After a few hundred seconds, 
the value of the angular velocity becomes quite insensitive to 
the value of K. The initial values of the amplitude of the 
mode and of the angular velocity of the star are, respectively, 
ao = 10 -6 and flo = Ok • 

the evolution, the angular velocity fl does not depend 
on the value of K, and, consequently, does not depend 
on the saturation value of a. As can be seen in Fig. |3 
already at t — to — 1000 s, for values of K ranging from 
to 10 s , the angular velocities are not much different; for 
t — t ~ 1 year any difference becomes negligibly small. 
This contrasts with the results obtained in Ref. Q , where 
the value of f2 after about one year depends critically on 
the choice of a sa t (see case K = — 2 below). 

After about one year of evolution, when the dissipative 
effect of viscosity becomes dominant and starts damp- 
ing the mode, the angular velocity of the star reaches 
values consistent with observational results. Indeed, as- 
suming that initially the star rotates with the maximum 
allowed angular velocity, f2o = — (2/3)\/7rGp, one 
obtains from Eq. l|3f)|) that f^one year — 0.042£lx, in good 
agreement with the inferred initial angular velocity of the 
fastest pulsars associated with supernova remnants. 

From Eqs. (gSJl and EHJ follows that, for t > few x 10 2 
seconds, the quantity (n/Qo)a 2 does not depend on time, 
just on K, namely, (tt/flo)a 2 ~ 8/(K+2), implying that 
the fraction of the initial angular momentum of the star 
that is transferred to the r-mode is just a function of K, 



5 {2) J/Jo 



if + 5/4 
K + 2 ' 



for t > few x 10 2 seconds, (39) 
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where we have taken Jo ~ JSlo ■ Thus, for K 1, most 
of the initial angular momentum of the star 7J7 is trans- 
ferred to the r-mode perturbation and, consequently, 
almost none is carried away by gravitational radiation 
emission. On the other hand, for K = —5/4, the angular 
momentum of the perturbation 5^ 2 'J remains zero dur- 
ing the evolution and, consequently, all the initial angu- 
lar momentum of the star IQq is available to be radiated 
away by gravitational radiation (see Figs. and |2J|. 



ferential rotation Att as the weighted variance of Q |F 



% 0.4 




3 4 5 

logi [(t-to)/sec] 

FIG. 3: Time evolution of the physical angular momentum 
of the r-mode perturbation J for different values of K. 
For K 3> 1, most of the initial angular momentum of the 
star Jo ~ Jf2o is transferred to the r-mode perturbation and, 
consequently, almost none is carried away by gravitational 
radiation emission. The initial values of the amplitude of the 
mode and of the angular velocity of the star are, respectively, 
do = 10~ 6 and fio = Qk- 
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FIG. 4: Time evolution of the total angular momentum of 
the star J for different values of K. For small values of K, a 
significant part of the initial angular momentum of the star is 
radiated away by gravitational radiation. The initial values 
of the amplitude of the mode and of the angular velocity of 
the star are, respectively, oiq = 10~ 6 and fio = Qk- 

During the nonlinear evolution, the fluid develops a 
strong differential rotation. Let us define the average dif- 



(Afi) 2 = 



Jpr 2 



?(<5( 2 V - n dr ) 2 dv 



J pr 2 sin 2 OdV 



(40) 



where the average angular velocity fldr, characterizing 
the drift of fluid elements along stellar latitudes, is given 
by 



5™ J _ J pr 2 sin 2 eS^v'f'dV 
I ~ J pr 2 sin 2 OdV 



(41) 



For the polytropic model we have been considering, 
Eqs. (|40[l and (|41|l yield for the average differential rota- 
tion the following expression: 



AO = \a 2 nQ 



15 , , > IH 

— (2AK 2 + 56K + 35) 
56 v J 2 



(4X + 5) 5 



1/2 



(42) 



where H = J Q p(r)r 8 dr / (M R 6 ) = 0.01. As can be seen 
from Fig. after a few hundred seconds the average dif- 
ferential rotation increases rapidly, saturating at high val- 
ues relatively to the initial angular velocity of the star. 
Note that the average differential rotation never vanishes 
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FIG. 5: Time evolution of the average differential rotation Afi 
for different values of K. After a few hundred seconds the 
average differential rotation increases rapidly, saturating at 
high values relatively to the initial angular velocity of the star. 
The initial value of the amplitude of the mode is «o = 10~ 6 . 

for any value of K. This is a consequence of the fact that 
the second-order velocity field (|10fl has two components, 
one induced by first-order quantities and another, fixed 
by initial data, which is a pure second-order effect. Since 
these two terms cannot cancel each other, a velocity drift 
of fluid elements along stellar latitudes is an unavoidable 
feature of the nonlinear r-mode pulsation. 

Let us now turn to the case K < —5/4. As already 
mentioned above, at a certain point of the evolution, the 
r-mode has to be saturated by hand in order to avoid 
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that the total angular momentum of the star becomes 
negative. Thus, during a first stage of the evolution, a 
and SI are determined from Eqs. 1)28(1 and Q29JI and during 
a second stage, a — a sat ^ \J — 3/[(AK + 5)Q] and f2 is 
determined by Eq. (|30|) . The fact that a saturation value 
for a has to be fixed by hand introduces an element of 
arbitrariness which did not exist in the case K —5/4. 
In particular, by choosing appropriately a sa t it is possible 
to arrange that the final value of f2 is consistent with 
observational data and, additionally, that a significant 
part of the initial angular momentum of the star is carried 
away by gravitational waves. 

For -2 < K < -5/4 and K < -2, the solution of the 
system of equations (l2~%|l and flty is given by Eqs. 
and (|36[1 . During the second stage of evolution, a = 
a sat ^ v/-3/[(4X + 5)g] and Q(t) is given by 



0.030aL 



1- 



l(4K + 5)Qa 2 sat 

6 -i -1/6 



(43) 



where is the time at which occurs the transition from 
the first to the second stage of evolution and t— t* is given 
in seconds. During the two stages of evolution we have 
just described, the total angular momentum of the star 
J(t) decreases due to gravitational- wave emission. At the 
end of the second stage of the evolution, it is given by 



J(t f 



inal ) 



•/(, 



i + i(4K 



5)QaL 



^{t final) 



(44) 



where we have taken Jo — I^o and S^(t) is given by 
Eq. Finally, let us point out that, for K < -2, 

during the first stage of evolution, the angular velocity 
of the star increases. 

The case K = — 2 has to be treated separately. It 
corresponds to the model studied in detail in Ref. p ■ The 
system of equations (|28|l and l|29|) yields the solution Vt = 
VLq and a = ao exp{— (t — io)/ T Gfl}- If the initial angular 
velocity is chosen to be flo — £Ik, then tqr = —37.1s, 
implying that the perturbation grows exponentially from 
the initial amplitude ao = 10 -6 to values of the order 
unity in just about 500 s |{J. After this short initial period 
in which f2 is constant and a grows exponentially, the 
amplitude a has to be forced, by hand, to take a certain 
saturation value a sa t ^ Q^ 1 ^ 2 — 3.26, and the angular 
velocity of the star is then given by Eq. (|43|l . with K = 
—2 and fl(t*) = f^o- The final angular momentum of 
the star is given by Eq. I|44|) with K = — 2. As can be 
seen from these equations, the final angular velocity and 
momentum of the star depend critically on the saturation 
value of the mode's amplitude a sa t] for instance, after one 
year of evolution, for a sa t = 1 and Qq = Qk one obtains 
{1 ~ O.lOif and J ~ 0.09J , while for a sat = 10~ 3 the 
angular velocity and momentum are, respectively, £1 ~ 
0.9fl K and J~ 0.9J . 



V. CONCLUSIONS 

In this paper we have studied the role of differential 
rotation in the evolution of the r-mode instability. We 
have adopted the simple model of Ref. , with two mod- 
ifications: (i) the physical angular momentum of the r- 
mode perturbation includes not only the canonical an- 
gular momentum but also a piece linear in second-order 
quantities, corresponding to differential rotation inducing 
large scale drifts of fluid elements along stellar latitudes; 
(ii) the evolution equations are deduced just from angu- 
lar momentum considerations. The first modification is 
a quite important one, resulting from the fact that dif- 
ferential rotation is an unavoidable kinematic feature of 
r- modes ^jj- The presence of this differential rotation 
implies that at second order in the mode's amplitude r- 
modes do not preserve vorticity of fluid elements. This in 
turn implies that the canonical angular momentum is not 
the full angular momentum at second order; one should 
also include a part linear in the second-order Lagrangian 
change in velocity, which is related to conservation of 
circulation in the fluid. The second modification, less 
significant, leads to a system of equations for the evolu- 
tion of the r-mode instability that differs from the ones 
of Ref. |6j just in a quantity of the order of Qa 2 . A de- 
tailed justification for this modification can be found in 

Ref. Eg. 

Within this model, we have derived a system of two 
first-order, coupled, differential equations (|2l)|) and QZ7\ . 
determining the time evolution of the amplitude of the 
r-mode a(t) and the angular velocity of the star Q(t). 
For the gravitational and viscous timescales appearing 
in these equations, we have used the expressions derived 
in Refs. for the simple case of r-modes of a newly 
born, hot, rapidly rotating neutron star with shear and 
bulk viscosity. In this case, the driving effect of the grav- 
itational radiation reaction overcomes the damping effect 
of shear and bulk viscosity for about one year, while the 
temperature of the star decreases from about 10 10 K to 
about 10 9 K. During this period of time, in which the r- 
mode instability is active, the gravitational timescale is 
much smaller than the viscous timescale, tgr. *C Ty , and, 
therefore, the evolution of a and £1 can be determined, in 
a good approximation, by the system of equations l|28|l 
and JUJl. 

The system of equations (I28f) and (|29l) contains a pa- 
rameter K, which is fixed by initial data and gives the 
initial amount of differential rotation associated with the 
r-mode. The specific case K — —2, for which the physi- 
cal angular momentum of the r-mode perturbation coin- 
cides with the canonical angular momentum, was studied 
in great detail in Ref. There is not, however, to our 
knowledge, any physical condition that forces K to take 
such a particular value K = —2. Therefore, we have 
solved, both numerically and analytically, the system of 
equations (gHJ and lO for arbitrary \K\ <C 10 13 . This 
upper limit for K results from the fact that one wishes 
to impose the condition that the initial absolute value 
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of the physical angular momentum of the perturbation 
J(to)\ is much smaller than the angular momentum 
of the unperturbed star I^Iq . In terms of the initial drift 
of fluid elements along stellar latitudes, this condition 
implies that (vdr)o at the equator of the star. It 

could be argued that some dissipative mechanism, such as 
bulk viscosity or magnetic coupling, would reduce differ- 
ential rotation in a newly born, hot, rapidly rotating neu- 
tron star, before the r-mode instability sets in, implying 
that K ~ 0. However, even in that case, a residual dif- 
ferential rotation of r-modes would be present in the star 
and the average differential rotation Af2 would increase 
exponentially as the amplitude of the r-mode grows, sat- 
urating at values relatively high as compared with the 
initial angular velocity of the star. Finally, in what con- 
cerns the admissible values of K, let us point out that for 
K < —5/4, as the amplitude of the mode grows due to 
the gravitational-radiation instability, the total angular 
momentum of the star decreases and eventually becomes 
negative. To avoid this unphysical situation, the growth 
of the mode's amplitude has to be stopped by hand at a 
saturation value a sa t ^ yJ—3/[(4:K + 5)Q] and integra- 
tion has then to be carried on with a new set of equations. 
The fact that a saturation value for a has to be fixed by 
hand introduces an element of arbitrariness into the so- 
lution, permitting, for instance, that agreement between 
the predicted final value of the angular velocity of the 
star and the value inferred from astronomical observa- 
tions can always be achieved by simply fine tuning the 
value of a sa t- One would wish, of course, that such an 
agreement, even qualitative, would arise in a natural way, 
without fine tuning the parameters of the model. For 
this reason, we have concentrated most of our attention 
on the case K ^ —5/4, in which such arbitrariness does 
not exist. 

From the exact analytical solution l|33|) and (|36[1 of the 
system of equations l|28|) and (|29|) one can extract several 
conclusions. 

First, the amplitude of the r-mode saturates in a natu- 
ral way a few hundred seconds after the mode instability 
sets in. The saturation amplitude depends on the param- 
eter K, namely, a sa t oc (K + 2)~ x / 2 . Therefore, if the 
initial differential rotation of r-modes is small (K ~ 0), 
then the r-mode saturates at values of the order of unity. 
On the other hand, if the initial differential rotation is 
significant (K 3> 1), then the saturation amplitude can 
be as small as 1CP 3 — 10~ 4 . These low values for the 
saturation amplitude of r-modes are of the same order 
of magnitude as the ones obtained in recent investiga- 
tions on wind-up of ma gne tic fields [9j and on nonlinear 
mode- mode interaction |l6j . 

Second, the value of the angular velocity of the star 



becomes, after a short period of evolution [t 3> few x 10 2 
s), very insensitive to the value of the parameter K, i.e., it 
becomes insensitive to the saturation value of the mode's 
amplitude. From Eq. H36|) for t = t(fi), it can be easily 
obtained that the angular velocity after about one year 
of evolution of the r-mode instability is 0.042f2x (for any 
K), in good agreement with the inferred initial angular 
velocity of the fastest pulsars associated with supernova 
remnants. 

Finally, the value of the physical angular momentum 
S^J of the r-mode perturbation tends, after a short pe- 
riod of evolution (t ^> few x 10 2 s), to a constant value 
given approximately by [(K + j)/(K + 2)] Jo. Thus, for 
K ^> 1 most of the initial angular momentum of the star 
is transferred to the r-mode and, consequently, almost 
none is carried away by gravitational radiation. That 
the r-mode absorbs most of the initial angular momen- 
tum of the star is explained by the fact that the mode 
develops a very strong differential rotation. Indeed, the 
average differential rotation, given by Eq. (|42[1 , increases 
exponentially, saturating after a few hundred seconds of 
evolution at very high values relatively to the initial an- 
gular velocity of the star. On the other hand, for values 
of K ~ —5/4, the transfer of angular momentum to the 
r-mode is less strong and, consequently, more angular 
momentum is available to be carried away by gravita- 
tional radiation. 

As we have seen, differential rotation introduces in the 
evolutionary picture of the r-mode instability new and 
somewhat unexpected features. This differential rotation 
has a kinematic origin, being induced by terms quadratic 
in the velocity field of the linear r-mode. However, differ- 
ential rotation of r-modes could also be induced by the 
gravitational-radiation reaction force. Indeed, within a 
toy model of a thin spherical shell of a rotating incom- 
pressible fluid it was shown that differential rotation of 
r-modes is also driven by a gravitational radiation force, 
leading to the conjecture that in real stars one could ob- 
serve a similar behavior j2jj . The existence of this differ- 
ential rotation induced by gravitational radiation and its 
influence on the evolution of the r-mode instability will 
be addressed in a future publication pi) . 
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